Add history outputs: Crop biomass and LAI per harvest#3633
Add history outputs: Crop biomass and LAI per harvest#3633samsrabin merged 28 commits intoESCOMP:b4b-devfrom
Conversation
Was using repr_grainc_to_food_patch + repr_grainc_to_seed_patch, but those don't get updated until after the harvest in CropPhenology, which is where the last CropPhaseTransitionBiomass() is called. They also, obviously, only include grain. Now using reproductivec_patch.
|
Adding "next" label until we get a reviewer. |
|
@slevis-lmwg can you take this one one, I think you have the most familiarity with the crop model. Is it helpful to go through the PR with @samsrabin ? |
slevis-lmwg
left a comment
There was a problem hiding this comment.
@samsrabin thank you for going over this work with me. I'm interested in looking again when you resolve the comment that you posted: https://github.com/ESCOMP/CTSM/pull/3633/changes#r2677039517
Also we agreed that this is fine in b4b-dev because it changes field-lists and not answers.
|
It's also possible for harvest to occur without having reached anthesis, or even emergence.
This will require moving |
|
Blocked by #3713. |
|
#3713 is merged, so this is unblocked. |
…leafout and grainfill.
slevis-lmwg
left a comment
There was a problem hiding this comment.
Thank you @samsrabin
Description of changes
Adds a bunch of history outputs with the
mxharvestsdimension (max number of crop harvests per year; 2).MAX_TLAI_PERHARVgives the maximum total LAI seen across the growing season.20 other outputs give the biomass for five C pools (leaf, grain, stem, fine root, coarse root) at three crop phase transitions (emergence, anthesis, maturity) and harvest (which may not happen at maturity).
Also adds an option to suppress a warning about low
gddmaturity, which I used during the CRU-JRA crop calendar work.Specific notes
Contributors other than yourself, if any: None
CTSM Issues Fixed (include github issue #):
Are answers expected to change (and if so in what way)? No
Any User Interface Changes (namelist or namelist defaults changes)? No
Does this create a need to change or add documentation? Did you do so? Will need to update non-FATES output list; no.
Testing performed
✅ New unit tests
Pass.
✅ aux_clm
Resolved with a7e66aa.aux_clmreveals what seem to be threading problems withMAX_TLAI_PERHARV, which failsCOMPARE_base_restinER[PS]_Ly3_P64x2tests but not anERS_Ly3test.✅ Sense checks
To do
Add HARVORGANC equivalent but JUST for grain. Will match what I have there for now, because grain is all CLM simulates by default. But it'll differ with the options underneath. (Do not name GRAIN unless you fix LPREPRSTRUCT test to handle this.)People should just use existingGRAINC_TO_(FOOD|SEED)_PERHARVoutputs.Add support forSee above.for_testing_use_second_grain_pooland/orfor_testing_use_repr_structure_pooltrue (or throw error in that case).